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1. Introduction 

In this paper we study the following SPDE with obstacle (in short OSPDE): 

' du t (x) = di (a i) j(t,x)d j u t (x)) dt + d i gi(t,x,u t (x),'Vut(x))dt + f(t,x,u t (x), Vu t (x))dt 
+00 

+ \ hj(t, x, u t (x), \7ut{x))dBl + v(t, dx), 
u t > S t , 

, ^0 = £ ■ 

(1) 
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where a is a time-dependant symmetric, uniformly elliptic, measurable matrix denned on 
some open domain O C M. d , with null Dirichlet condition. The initial condition is given 
as no = £, a L 2 (0)— valued random variable, and f,g = (gi, ■■■,gd) an d h = (h±, ...hi, ...) 
are non-linear random functions. Given an obstacle S : £1 x [0, T] x O — > R, we study the 
obstacle problem for the SPDE (JXJ) , i.e. we want to find a solution of (JXJ) which satisfies 
"u > S" where the obstacle S is regular in some sense and controlled by the solution of a 
SPDE. 

In recent work |10| we have proved in the homogeneous case, existence and uniqueness of 
the solution of equation (pQ) with Dirichlet boundary condition under standard Lipschitz 
hypotheses and L 2 -type integrability conditions on the coefficients. Moreover in still in 
the homogeneous case, we have obtained a maximum principle for local solutions. In these 
papers we have assumed that a does not depend on time and so many proofs are based on 
the notion of semigroup associated to the second order operator and on the regularizing 
property of the semigroup. The aim of this paper is to extend all the results to the non 
homogeneous case. 

Let us recall that the solution is a couple (u,u), where u is a process with values in the 
first order Sobolev space and v is a random regular measure forcing u to stay above S and 
satisfying a minimal Skohorod condition. In order to give a rigorous meaning to the notion 
of solution, inspired by the works of M. Pierre in the deterministic case (see \2?>\ 124]). we 
introduce the notion of parabolic capacity. We construct a solution which admits a quasi 
continuous version hence defined outside a polar set and use the fact that regular measures 
which in general are not absolutely continuous w.r.t. the Lebesgue measure, do not charge 
polar sets. 

There is a huge literature on parabolic SPDE's without obstacle. The study of the LP — norms 
w.r.t. the randomness of the space-time uniform norm on the trajectories of a stochastic 
PDE was started by N. V. Krylov in |16| (see also Kim |14|). for a more complete overview 
of existing works on this subject see [HJ [9] and the references therein. Let us also men- 
tion that some maximum principle have been established by N. V. Krylov |17| for linear 
parabolic spde's on Lipschitz domain. Concerning the obstacle problem, there are two ap- 
proaches, a probabilistic one (see |20[ I15| ) based on the Feynmann-Kac's formula via the 
backward doubly stochastic differential equations and the analytical one (see |12 [ I22 [ [2~T]) 
based on the Green function. 

The main results of this paper are first an existence and uniqueness Theorem for the 
solution with null Dirichlet condition and a maximum principle for local solutions. This 
yields for example: 

Theorem 1. Let (Mf)t>o be an Ltd process satisfying some integrability conditions, p > 2 
and u be a local weak solution of the obstacle problem (Op. Assume that dO is Lipschitz 
and u < M on dO, then for all t £ [0, T] : 

E\\(u- M) + ||^ )0o;t < k (p, t) C(S, /, g, h, M) 

where C{S, f,g,h,M) depends only on the barrier S, the initial condition £, coefficients 
f, g, h, the boundary condition M and k is a function which only depends on p and t, 
II ■ ||oc oo-t i> s the uniform norm on [0,t] x O. 
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2. Hypotheses and preliminaries 
2.1. Settings 

Let O be an open bounded domain in W 1 . The space L 2 (O) is the basic Hilbert space of 
our framework and we employ the usual notation for its scalar product and its norm, 

i 

2 



(u,v) = / u(x)v(x)dx, \\u\\ = I / u (x) dx 
Jo \Jo 

In general, we shall extend the notation 

(u, v) = J u(x)v(x)dx, 
Jo 

where u, v are measurable functions denned on O such that uv G L l (0). 

The first order Sobolev space of functions vanishing at the boundary will be denoted as 

usual by Hq (0) . Its natural scalar product and norm are 



i 

|2\ 2 



d „ 

(u, w)/fi(o) = (u, v) + ^ / diu (x) div (x) dx, = M|^i|)§ + || Vt 

i=i J° 

We shall denote by i?/ oc (0) the space of functions which are locally square integrable in 
O and which admit first order derivatives that are also locally square integrable. 
Another Hilbert space that we use is the second order Sobolev space Hq(0) of functions 
vanishing at the boundary and twice differentiable in the weak sense. 

We consider a sequence ((S*(t))t>o)jgN* of independent Brownian motions defined on a 
standard filtered probability space (0, J 7 , (J r j)t>o,-P) satisfying the usual conditions. 
Let a be a measurable and d x d symmetric matrix defined on M + x O. We assume that there 
exist positive constants A, A and M such that for all r\ G W 1 and almost all (i, x) S ]R + x O: 

Mv\ 2 < a ij{t, x)rfrf < A\r]\ 2 and \aij(t,x)\ < M. (2) 

Let A = {(t, x, s, y) G E+ x O x E+ x O; t> s}. We denote by G : A E+ the weak 
fundamental solution of the problem 

d 

d t G(t,x;s,y) - } djQi,j(t, x)djG(t, x; s, y) =0 (3) 

i=l 

with Dirichlet boundary condition G(t,x; s,y) = 0, for all (t,x) G (s,+oo) x dO. 

We consider the quasilinear stochastic partial differential equation (1) with initial condition 

u(0, •) = £(•) and Dirichlet boundary condition u(t,x) = 0, V (t, x) G M + x dO. 

We assume that we have predictable random functions 

JiR+xfixOxlxR^l, 

g = (g 1 ,...,g d ) :R + x!]xOxlxK d ->K d , 

h = (h u ..., hi,...) :M+ xflx O xRxR d -> M N *, 
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In the sequel, | ■ | will always denote the underlying Euclidean or I 2 — norm. For example 



+00 

\h(t,uj,x,y,z)\ 2 = ^ \hi(t,uj,x,y,z)\ 2 . 
i=i 

Assumption (H): There exist non- negative constants C, a, f3 such that for almost all 
uj, the following inequalities hold for all (x,y,z,t) £ O x R x R d x R + : 

I. \f(t,u),x,y,z) - f(t,u,x,y',z')\ < C{\y - y'\ + \z- z'\), 

2 - (EiLi \9i{t,u,x,y,z) - gi(t,uj,x,y',z')\ 2 )2 < C\y - y'\ + a\z - z'\, 

3. (\h(t,u,x,y,z) -h^u.x.y'.z')^ < C\y - y'\ + P\z - z'\, 

4. the contraction property: 2a + f3 2 < 2A. 

Remark 1. This last contraction property ensures existence and uniqueness for the solution 
of the SPDE without obstacle (see 

Moreover for simplicity, we fix a terminal time T > 0, we assume that: 
Assumption (I): 

£ E L 2 (il x O) is an J-q — measurable random variable 

/(.,-, -,0,0) :=f°€L 2 ([0,T]xnxO;R) 

g(-, ■, ■, 0, 0) := g° = (<??, ...,g° d ) € L 2 ([0, TjxQxO; R d ) 

h(-,;;0,0) : = h° = (h?, h?, ...) G L 2 ([0,T] x x 0;M N *). 

We denote by TLt the space of Hq(0)— valued predictable L 2 (0)— continuous processes 
( u t)te[o,T] which satisfy 

\\u\\t = E sup |j ut || 2 +-E 1 / || Vut || 2 (it < +00. 
te[o,T] 

It is the natural space for solutions. 

The space of test functions is denote by T> = C^°(R + ) (g>C 2 (0), where C£°(K + ) is the space 
of all real valued infinitely differentiable functions with compact support in R + and C 2 {0) 
the set of C 2 -functions with compact support in O. 

Main example of stochastic noise 

Let W be a noise white in time and colored in space, defined on a standard filtered prob- 
ability space (p., J-, (J-t)t>o, P) whose covariance function is given by: 

Vs, teR+, Vx,y£(D, E[W(x, s)W(y, t)] = S(t - s)k(x, y), 

where fc:Ox04 R + is a symmetric and measurable function. 
Consider the following SPDE driven by W: 

d d 

du t (x) = ( ^2 diaij(t,x)djUt(x) + f(t,x,ut(x),Vut(x)) + 22^i9i(t,x,Ut(x),Vu t (x)))dt 

i,j=l i=l (4) 

+ h(t, x, u t (x),X7u t (x)) W(dt, x), 
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where / and g are as above and h is a random real valued function. 

We assume that the covariance function k defines a trace class operator denoted by K in 
L 2 (0). It is well known that there exists an orthogonal basis (ej)j g N* of L 2 (0) consisting 
of eigenfunctions of K with corresponding eigenvalues (Aj)j g N* such that 

+00 

Ai < +00, 

i=l 

and 

+00 

1=1 

It is also well known that there exists a sequence ((.B l (i))t>o)ieN* of independent standard 
Brownian motions such that 

+00 

,1/2, 



W(dt,-) = ^A i 1/2 e lJ B i (cft). 



i=l 

So that equation Q is equivalent to equation ([TJ without obstacle and with h = (/ii)igN* 
where 

Vi e N*, hi(s,x,y,z) = \/\h(s,x,y,z)ei(x). 
Assume as in |26] that for all i E N*, He^loo < +00 and 



i=l 

Since 

1 /+00 



(^[/i(i,o;,x,y,z) - /i(t,w,x,y',z')| 2 ) 2 < ( X^H^H^ J x ' y ' z ) ~ ^(*, x, y', z') 



a=l 



2 



/i satisfies the Lipschitz hypothesis (H)-(ii) if h satisfies a similar Lipschitz hypothesis. 



2.2. Parabolic potential analysis 

In this section we will recall some important definitions and results concerning the obstacle 
problem for parabolic PDE in |23| and |24| . 

K. denotes L°°([0, T];L 2 (0)) n L 2 ([0, T];Hq(0)) equipped with the norm: 

II II 2 II II 2 1 11 ||2 

II v \\K — II v llL°°([0,T];L 2 (O) + II V Wl 2 ([0,T};H^(O)) 

> II vt II 2 + / 

te[o; 

C denotes the space of continuous functions on compact support in [O,T[x0 and finally: 



sup II v t || 2 + / (II v t || 2 + II Vv t || 2 ) dt. 

'efo,T[ Jo 



W = We L\%T]-Hl(0))- ^ G L 2 ([0,T];H~\O))}, 



6 



L. Denis, A. Matoussi and J. Zhang 



II 9 ll 9 ll 2 

endowed with the norm|| <p \\fo=\\ <p \\ L 2 [[0)T] . H i io)) + II ^ \\V{[o,T];H-i(0)y 
It is known (see |18| ) that W is continuously embedded in C([0, T]; L (0)), the set of 
L 2 (0)-valued continuous functions on [0, T]. So without ambiguity, we will also consider 
W T = W G W; = 0}, W + = y € W; 9? > 0}, W + = W T n W+ 

We now introduce the notion of parabolic potentials and regular measures which permit 
to define the parabolic capacity. 

Definition 1. An element v £ K, is said to be a parabolic potential if it satisfies: 



V<^ E W, 



T > 



Jo ~^Et ,Vt)dt + Jo £ ( ( Pt,v t )dt>0. 



We denote by V the set of all parabolic potentials. 
The next representation property is crucial: 

Proposition 1. (Proposition 1.1 in \21$ ) Let v 6 V, then there exists a unique positive 
Radon measure on [0,T[xO, denoted by v v , such that: 



v^€W T nc, (_^ |Vt )dt+ s(ip t ,v t )dt= i ' ip(t, x )dv v . 



rT r T 
at Jo JO JO 

Moreover, v admits a right- continuous (resp. left- continuous) version v {resp. v) : [0, T] \— > 
L 2 {0) . 

Such a Radon measure, u v is called a regular measure and we write: 

v v = — + Av. 
at 

Remark 2. As a consequence, we can also define for all v € V: 

v T = \imvt € L 2 (0). 

tyr 

Definition 2. Let K C [0,T[xO be compact, v € V is said to be v— superior than 1 on 
K, if there exists a sequence v n 6 V with v n > 1 a.e. on a neighborhood of K converging 
to v in L 2 ([0,T];H^(O)). 

We denote: 

'5^K = {v € V; v is v — superior to 1 on K}. 

Proposition 2. (Proposition 2.1 in [24J) Let K C [0, T[xO compact, then S"k admits a 
smallest Vr € V and the measure v v K whose support is in K satisfies 



[ f dv v K =iTd{l [ dv v ;vey K }- 
Jo Jo v ^ v Jo Jo 



[ dvl = inf { f 
<o v&v Jo JO 

Definition 3. (Parabolic Capacity) 



Let K C [0, T[xO be compact, we define cap(K) = J^jdv 1 ^; 

let O C [0,T[xO be open, we define cap(0) = sup{cap(K); K C O compact); 

for any borelian E C [0,T[xO, we define cap(E) = inf {cap(0); O D E open}. 
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Definition 4. A property is said to hold quasi- everywhere (in short q.e.) if it holds outside 
a set of null capacity. 

Definition5. (Quasi- continuous) 

A function u : [0, T[xO — > R is called quasi- continuous, if there exists a decreasing sequence 
of open subsets O n of [0,T[xO with: 

1. for all n, the restriction of u n to the complement of O n is continuous; 

2. lim n _> +00 cap (O n ) = 0. 

We say that u admits a quasi- continuous version, if there exists u, quasi- continuous such 
that u = u a.e.. 



The next proposition, whose proof may be found in [23j or [24j shall play an important 
role in the sequel: 

Proposition 3. Let K C O a compact set, then Vi € [0, T[ 

cap({t} xK) = X d (K), 
where Xd is the Lebesgue measure on O. 

As a consequence, if u : [0,T[xO — > R is a map defined quasi- everywhere then it defines 
uniquely a map from [0,T[ into L 2 (0). In other words, for any t € [0,T[, Ut is defined 
without any ambiguity as an element in L 2 (0). Moreover, if ' u € V, it admits version u 
which is left continuous on [0, T] with values in L 2 (0) so that ut = Ur- is also defined 
without ambiguity. 

Remark 3. The previous proposition applies if for example u is quasi- continuous. 

Proposition 4. (Theorem III.l in \21$) If(f€ W, then it admits a unique quasi- continuous 
version that we denote by (p. Moreover, for all v 6 V, the following relation holds: 

rT 

(pdv v = I (-dt<p,v)+£((p,v)dt + (<pT,vr)- 



[0,T[xO Jo 

We end this section by a convergence lemma which plays an important role in our approach 
(Lemma 3.8 in [23]): 

Lemma 1. Ifv n E V is a bounded sequence in K and converges weakly to v in L 2 ([0, T]; Hq(O)); 
if u is a quasi- continuous function and \u\ is bounded by a element in V. Then 

lim [ T [ udv vn = [ T [ udu v . 
/o Jo Jo Jo 



Remark 4. For the more general case one can see WM Lemma 3. 



3. Quasi-continuity of the solution of SPDE without obstacle 

We consider the SPDE without obstacle: 

du t (x) = di (a itj (t,x)dju t (x) + gi(t,x,u t (x),X7u t (x))) dt + f(t,x,u t (x),Vu t (x))dt 

+oo 

+ hj (t, x, u t (x) , Vu t (x))dB{ , (5) 

3=1 
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As a consequence of well-known results (see for example [9j, Theorem 11), we know that 
under assumptions (H) and (I) , SPDE ([5]) with zero Dirichlet boundary condition, admits 
a unique solution in ~Ht, we denote it by f,g,h), moreover it satisfies the following 
estimate: 

E[\\u\\l] < cE 



Uf + l (llA°H 2 + Hlrf + IH* 



0|||2 



dt 



(6) 



The main theorem of this section is the following: 



Theorem 2. Under assumptions (H) and (I), u = U(£, f,g,h) the solution of SPDE (0j 
admits a quasi- continuous version denoted by u i.e. u = u dP®dt®dx—a.e. and for almost 
all w € Q,, (t,x) — > Ut(w,x) is quasi- continuous. 

Before giving the proof of this theorem, we need the following lemmas. The first one is 
proved in |24j . Lemma 3.3: 

Lemma 2. There exists C > such that, for all open set $ C [0,T[xO and v £ V with 
v > 1 a.e. on 



cap'd < C || v ||£ 



Let k = k(u, u + (0)) be defined as following 

k = ess'mf{v 6?; v > u a.e. and v(0) > u + (0)}. 

One has to note that k is a random function. From now on, we always take for k the 
following measurable version 

ft = supv n , 

n 

where (v n ) is the non-decreasing sequence of random functions given by 



dt t \ t t, ^ 



Lv? + n(v? - utY 



< = n + (0). 

From F.Mignot and J.P.Puel |21| . we know that for almost all w G fi, v n (w) converges 
weakly to v(w) = k(u(w), u + (0){w)) in L 2 ([0, T];Hq(0)) and that v > u. 

Lemma 3. We have the following estimate: 

E\\k\\1<C (e (I u+ || 2 +E || n || 2 +E C \\ f t || 2 + || \g*\ f + || \h\\ f dt) , 



o 



where C is a constant depending only on the structure constants of the equation. 



Thanks to (|2|), the proof of Lemma 3 in |10| can be easily extended to the case of non- 
homogeneous operator. 

Proof of Theorem [2} First of all, we remark that we only need to prove this result in 
the linear case, namely we consider that /, g and h only depend on t, x and uj. Then, we 
approximate the coefficients, the domain and the second order operator in the following 
way: 







1. We mollify coefficients a^j and so consider sequences (a"j) n of C°° functions such 
that for all n G N*, the matrix a n satisfies the same ellipticity and boundedness 
assumptions as a and 

Vf <i,j < d, lim a",- = 7 - a.e. 

2. We approximate by an increasing sequence of smooth domains {O n ) n >\. 

3. We consider a sequence (£ n ) in C%°(0) which converges to £ in L 2 (0) and such that 
for all n, supp£ n C O n . 

4. For each % G N*, we construct a sequence of predictable functions (hf) in 

(L 2 (n) <g> C c ([0, +oo)) <8) C c °°(0)) which converges in L 2 oc (M+; ^ 2 (^ x C)) to h i such 
that for all n, suppw™ C O n and 

Vt>0, E[f \\hl s \\ 2 ds] < E[ [ \\h itS \\ 2 ds], 
Jo ' Jo 

so that ^ 

E[l \\\h n s \\\ 2 ds\<E[f \\\h s \\\ 2 ds] < +oo. 
Jo Jo 

5. We consider a sequence of predictable functions (/ n ) in (L 2 (il) (g) C c ([0, +oo)) <g> C%°(0)) 
which converges in Lf oc (R + ; L 2 (U x O)) to / and such that for all n, supp/ n C O n . 

6. Finally, let (g n ) be a sequence in (L 2 (fi) ® C c ([0, +oo)) ® C~(0)) which converges 
in L 2 oc (M+;L 2 (Sl x to #and such that for all n, supp# n C O n . 

For all n G N*, we put A n = {(t,x,s,y) G R+ x O n x R + x O n ; t > s}. We denote by 
G n : A n h-> M + the weak fundamental solution of the problem ([3]) associated to a n and 
O n : 

d 

d t G n (t,x;s,y) - ^ d i al j (t,x)d j G n {t,x-s,y) = Q (8) 
»,j'=i 

with Dirichlet boundary condition G n (t,x,s,y) = 0, for all (t, x) G (s, +oo) x <9C? n . 
In a natural way we extend G n on A by setting: G n = on A \ A" . 
We define the process u n by setting for all (t, x) G 1R+ x O: 



u?(.)= / G(t,-,0,y)e n (y)dy+ / / ■, s, y)^(y)dyds 
Jo Jo Jo 

+ ^2 / G ( t r,s,y)di, y g?(y)dyds ^ 

+oo ,. t ,. 

+ J2 G(t,-,s,y)hl s (y)dBl. 
~r, Jo Jo 



The main point is that there exists a subsequence of (G n ) n >i which converges everywhere 
to G on A, where G still denotes the fundamental solution of ([3]), see Lemma 7 in [9]. 
From Proposition 6 in [9], we know that u n G Ht is the unique weak solution of (|5j). G is 
uniformly continuous in space-time variables on any compact away from the diagonal in 
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time ( see Theorem 6 in PQ) and satisfies Gaussian estimates (see Aronson [2]), this ensures 
that for all n € N*, u 11 is P-almost surely continuous in (t, x). 

Moreover, since sequences (f n ) n , {g n )n an d {h n )n are uniformly bounded in L 2 -spaces, as 
a consequence of estimate ([6]), (u n ) n is bounded in T~Lt hence in L 2 ([0,T] x Q;Hq(0)), so 
that we can extract a subsequence (u nk )k which converges weakly in L 2 ([0, T] x Q; Hq(0)) 
and such that a sequence of convex combinations (u n ) of the form 

N n 

fc=i 

converges strongly to u in L 2 (tt x [0, T\; Hq (O)). It is clear that for all n € N*, (u n ) n is 

P— almost surely continuous in (t,x). 

We consider a sequence of random open sets 



n=p 



Let K n = n{±(u n+1 - u"), ±{u n+1 - u n ) + {0)) + K(-±(u n+l - u n ), ±(u n+l - u n )~(0)), 



from the definition of k and the relation (see 

k(M) < k(v,i; + (0)) +k(-i;,u-(0)), 

we know that K n satisfy the conditions of Lemma [2j i.e. K n € V et K n > 1 a.e. on thus 
we get the following relation 



cap(Q p ) < ^ cap (■!?„) < ^ || K n 



1 2 

Ik 



Thus, remarking that u n+1 - u n = U(£ n+1 - £ n , f n+1 - f n , g n+1 - g n , h n+1 - h n ), we apply 
LemmaOto ft(^(u" +1 -u n ), ±(u n+1 -u n )+(0)) and K (-^(u n+1 -«"), ±(u n+1 -u n )~{0)) 
and obtain: 



+oo ^ „T 

£[cap (G p )] < II Kn HI < 20^^(2 H n+l -e\\ 2 +E / || /f +1 - if 1 

+ iibr +1 -5 t n iii 2 + ni^ +1 -^iii 2 ^)- 



Then, by extracting a subsequence, we can consider that 

e ii e n+1 - e ii 2 +^ £ ii /r 1 - /r n 2 + n - <?ri n 2 + n - w n 2 * < ^. 



Then we take €„ = ^ to get 



±2^ 2C7n 4 
E[cap (6 P )] < ^ — -. 



2 r 
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Therefore 

lim E[cap (©„)] = 0. 

p— t+oo 

For almost all u) € f2, u u (uj) is continuous in (i, x) on (0p(u;)) c and ({t n (ci;)) n converges 
uniformly to u on (Q p (w)) c for all p, hence, u(oj) is continuous in (i, x) on (0 p (u;)) c , 
then from the definition of quasi-continuous, we know that u(uS) admits a quasi-continuous 
version since cap (0 p ) tends to almost surely as p tends to +oo. □ 



4. Existence and uniqueness result 

From now on, similarly to the homogeneous case studied in |10j . we make the following 
assumptions on the obstacle: 

Assumption (O): The obstacle S is assumed to be an adapted process, quasi-continuous, 
such that So < (, P-almost surely and controlled by the solution of a SPDE, i.e. Vt G [0, T], 

s t < S' t (10) 

where 5' is the solution of a linear SPDE 

dS' t (x) = ^(^•(t,x)5 i 5Kx)+^ t (x))dt + //(x)ctt + + EjS^, t ^)d^ n , 
S'(0) = S> [ > 

with S' € L 2 {VL x O) J'o-measurable, /', g' and h! adapted processes respectively in 
L 2 ([0,T] x n x £>;R), L 2 ([0,T] x n x C;M d ) and L 2 ([0,T] x x C;M N *). 

Remark 5. Here again, we know that S' uniquely exists and satisfies the following esti- 
mate: 



E sup || 5; || 2 +E f T || VS[ \\ 2 dt< CE \\\ S' || 2 + f (|| /; || 2 + || \g' t \ || 2 + || \ti t \ f)dt 
te[o:r] Jo I Jo 

(12) 

Moreover, from Theorem^ S' admits a quasi- continuous version. 

Let us also remark that even if this assumption seems restrictive since S ! is driven by the 
same operator and Brownian motions as u, it encompasses a large class of examples. 

We now are able to define rigorously the notion of solution to the problem with obstacle: 

Definition 6. A pair (u, v) is said to be a solution of the obstacle problem for |5[) if 

1. u € Ht and u(t, x) > S(t, x), dP ® dt ® dx — a.e. and uo(x) = £, dP <g> dx — a.e.; 

2. v is a random regular measure defined on [0, T) x O; 

3. the following relation holds almost surely, for all t € [0, T] and V(p € T>, 

pt ^ rt p 

(ut,tpt) - (€,<Po) - {u s ,d s ip s )ds - V] / / a i) j{s 1 x)d i u s (x)djLp s (x)dxds 
Jo i~Hi Jo Jo 

rt d rt 

= / (f s (u s , Vu s ),(p s )ds + y2 / (gl(u s ,Vu s ),ditp s )ds 
Jo l=1 Jo 

+00 „t /•' 

+ V / (hP a (u 8 ,Vu 8 ),(p 8 )dBi + / 
■ =1 Jo Jo 



-T 



00 rt 

ip s (x)u(dx,ds). (13) 



o 
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4- u admits a quasi- continuous version, u, and we have 



(u(s, x) — S(s, x))v(dx, ds) = a.s. 

'o Jo 

The first important result of this paper is: 

Theorem 3. Under assumptions (H) 7 (I) and (O), there exists a unique weak solution of 
the obstacle problem for the SPDE ([5p associated to (£, /, g, h, S). 

We denote by TZ(^, f, g,h, S) the solution of SPDE ^ with obstacle when it exists and is 
unique. 

Proof. As we have Ito's formula and comparison theorem for the solution of non homoge- 
neous SPDE ([5]), see Proposition 9 and Theorem 16 in [9], we can make the same proof as 
in the homogeneous case (see [10]). More precisely, we first establish the result in the linear 
case by following Section 5.2 in [10J. Then, we prove an Ito formula for the difference of 
two (linear) solutions of SPDE's with obstacle similarly to Section 5.4 in |10| and finally 
conclude thanks to a Picard iteration procedure as in Section 5.5 in |10] , □ 



We can also establish the following Ito formula and comparison theorem for the solution 
of SPDE ([5]) with obstacle. Here again, the proofs are the same as in |10] , 

Theorem 4. Let u be the solution of OSPDE (EJj and <J> : R + x R -)■ R be a function of 
class C ' . We denote by $' and the derivatives of $ with respect to the space variables 
and by ^ the partial derivative with respect to time. We assume that these derivatives are 
bounded and 0) = for all t > 0. Then for all t G [0, T], 

&(t,ut(x))dx + / / a i j(s,x)&'(s,u s (x))diU s (x)djU s (x)dxds = / &(0,£(x))dx 
o Jo Jo Jo 

+ [ [ ^-(s,u s (x))dxds + / (&(s,u s ), f s )ds - / / &'(s,u s (x))diU s (x)gi(x)dxds 
Jo Jo ds J f^Jo Jo 

+ V / (&(s,u s ),h j )dBi + -J2 / / $"(s,u s (x))(h jjS (x)) 2 dxds 
~iJo z ~[Jo Jo 



+ [ [ &(s,u s (x))v(dxds) , P - a.s. 

Jo Jo 



This Ito formula naturally leads to a comparison theorem, the proof being the same as 
in the homogeneous case (see Theorem 8 in |10]). More precisely, consider (u = 
^•(C 1 ) Z 1 ) 9 j h, S l ) the solution of the SPDE with obstacle 

du\{x) = di(aij(t, x)dju\{x) + gi{t, x, u\{x), Vu\{x))dt + / 1 (t, x, u\{x), \7u\{x))dt 

+oo 

+ ^2hj(t,x,uj(x),Vuj(x))dB J t + v l (x,dt) 
where we assume (£ , f 1 ,g,h) satisfy hypotheses (H), (I) and (O). 

We consider another coefficients / 2 which satisfies the same assumptions as / , another 
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obstacle S 2 which satisfies (O) and another initial condition £ 2 belonging to L 2 (S7 x O) 
and Fq adapted such that £ 2 > S 2 . We denote by (u 2 , v 2 ) = 1Z{^ 2 , f 2 ,g, h, S 2 ). 

Theorem 5. Assume that the following conditions hold 

1. C 1 <C 2 , dx®dP-a.e. 

2. f l (v}, Vu 1 ) < f(u l , Vu 1 ), dt ® dx <g> dP - a.e. 

3. S 1 < S 2 , dt®dx®dP - a.e. 

Then for almost all uj € Q, u l (t,x) < u 2 (t,x) q.e. 

5. Maximum principle for local solutions of the OSPDE 
5.1. L p,q -spaces 

For each t > and for all real numbers p, q > 1, we denote by L p ' q ([0, t] x O) the space of 
(classes of) measurable functions u : [0, t] X O — > R such that 



iff \l/P N 

(y |ti(s, :r)| p fix J ds 



1/9 



is finite. The limiting cases with p or q taking the value oo are also considered with the 
use of the essential sup norm. 

The space of measurable functions u : R_|_ — > L 2 (O) such that ||n|| 2 2 -t < 00 > f° r eacn £ > 0, 
is denoted by L 2 oc (R+; L 2 (C)) , where K+ denotes the set of non-negative real numbers. 
Similarly, the space Lf oc (R+; Hq (0)) consists of all measurable functions u : IR+ — > Hq (O) 
such that 

IMI 2 ,2;t + H Vu ll2,2;t < °°. 

for any t > 0. 

We recall that the Sobolev inequality states that 

IMI2* < cs || Vu|| 2 , 

for each u £ Hq (O) , where cs > is a constant that depends on the dimension and 
2* = if d > 2, while 2* may be any number in ]2, oof if d = 2 and 2* = 00 if d = 1. 
Finally, we introduce the following norm which is obtained by interpolation in L p,<? -spaces: 

IMI#;t = ll U ll2,oo;i V H U ll2*,2;t ' 

and we denote by L# ; t the set of functions u such that ||u||^. t is finite. Its dual space is a 
functional space: L*^. t equipped with the norm || ||^., t and we have 

t 

u (s, x) v (s, x) dxds < ||it||#. t IMI#-t , (14) 



'0 Jo 

for any u 6 £# ; t and v £ 



14 



L. Denis, A. Matoussi and J. Zhang 



5.2. Local solutions 

We define Hi oc = %/ c(C) to be the set of i^ oc (0)-valued predictable processes defined on 
[0, T] such that for any compact subset K in O: 

(E sup / u s (x) 2 dx + E I I | Vu s {x)\ 2 dxds \ < oo. 
o<s<tJk Jo Jk J 

Definition 7. We say that a Radon measure v on [0,T[xO is a local regular measure if 
for any non-negative <j) in C^°((D), (f)v is a regular measure. 

In (see Proposition 2.10), we have proved: 

Proposition 5. Local regular measures do not charge polar sets (i.e. sets of capacity 0). 
We can now define the notion of local solution: 

Definition 8. A pair (u, v) is said to be a local solution of the problem ([7]) if 

1. u € Hioc, u(t, x) > S(t, x), dP (g> dt <g> dx — a.e. and uq{x) = £, dP <8> dx — a.e.; 

2. v is a local random regular measure defined on [0, T[xO; 

3. the following relation holds almost surely, for all t 6 [0, T] and all ip € D, 



=(£>¥?o)+ / (u s ,d s ip s )ds - V] / / a i j{s 1 x)d i u s (x)dj(p s (x)dxds 
Jo jj Jo Jo 



d pt ft 

(gl(u s ,Vu s ),di<p s )ds + / (f s (u s ,Vu s ),(p s )ds 
l=1 Jo Jo 

+ 5^ / (h j a ( u s,'Vu 8 ),(p s )dB 3 8 + / / f s (x)u(dx,ds). 
~{Jo Jo Jo 

4- u admits a quasi- continuous version, u, and we have 

(u(s, x) — S(s, x))v(dx, ds) = 0, P — a.s. 



(15) 



'o Jo 

We denote by 1Z-i oc (£, f, 5S h, S) the set of all the local solutions (u, v). 



5.3. Hypotheses 

In order to get some L p — estimates for the uniform norm of the positive part of the solution 
of dl|), we need stronger integrability conditions on the coefficients and the initial condition. 
To this end, we consider the following assumptions: for p > 2: 
Assumption (HI2p) 

^ (iieii^ Il^°ll2, 2; ^ lli^°illl 2; ^ lli^°illl 2 ^) < ~- 

Assumption (OL): The obstacle S : [0, T] x Q, x O — > R is an adapted random field, 
almost surely quasi-continuous, such that Sq < £ P-almost surely and controlled by a local 
solution of an SPDE, i.e. Vi G [0,T], 

S t < S' t , dP ®dt®dx- a.e. 
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where S' is a local solution (for the definition of local solution see for example Definition 
1 in 0) of the linear SPDE 

f dS' t = LSldt + fldt + Y,Lidi9i,tdt + Et=iti jit dBi 
\ S'(0) = S' . 

Assumption (HIL) 

e\ \ax)\ 2 dx + E [ T [ (\f°(x)\ 2 + \g s (x)\ 2 + \h° s (x)\ 2 )dxds<oc, 

JK JO JK 

for any compact set K C O. 
Assumption (HOL) 



e! \S' \ 2 dx + E [ T [ (\fl(x)\ 2 + \g' t (x)\ 2 + \h' t (x)\ 2 )dxdt< 

JK JO JK 



DC 



for any compact set K C O. 
Assumption (HOoop) 

S> e x O) and E ((U/'L^ + (|| l</| 2 |L i00;T ) p/2 + (|| !^| 2 |L i0O;T ) p/2 ) < oo. 

As our approach is based on some estimates of u — S' that we obtain thanks to the Ito 
formula, we need to introduce the following functions: 

f(t, u), x, y, z) = f(t, u,x,y + S' t , z + VS' t ) - f'(t, u, x) 

g(t, u, x, y, z) = g(t, uj,x,y + S' t , z + VS' t ) - g'{t, u, x) 
h(t, u, x, y, z) = h(t, co, x, y + S' t , z + S/S' t ) — h'(t, lj, x). 

And we consider: 
Assumption (HD#p) 

^(||/i;t) p +(III5 o !i;t) p/2 +(Ii^ o ii;t) p/2 )<-- 

This assumption is fulfilled in the following case: 

Example 1. // || VS'\\l. T , \\f°\\* e . T , \\9°\\*q-t an d 1 1 1 1 S-T belong to L P (Q,P), and as- 
sumptions (H) and (HOoop) hold, then: 

f satisfies the Lipschitz condition with the same Lipschitz coefficients: 

\f(t,u,x,y,z) - f(t,uj,x,y' ,z')\ = \f(t,uj,x,y + S' t (x), z + VS' t (x)) + f'(t,u,x) 

- f(t, u, x, y' + S' t (x), z' + V^(x)) - f'(t, u, x) \ 

< C\y-y'\+C\z- z'\. 

f satisfies the integrability condition: 

\\f°\\l;T = ll/( 5, ' v ' s, )-/ / !i; T <!i/(5 / ,v5')ii; T +ii/ , ii; T 

- l|/°IL;T + C 'll' S ' / |le;T + C 'll V5 ' / |L;r+ l|/1L,oo;T ' 

And the same for g and h, which proves that (HD#p) holds. 
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5-4- The main results 

We now introduce the lateral boundary condition that we consider: 

Definition 9. If u belongs to T-Li oc , we say that u is non-negative on the boundary of O if 
u + belongs to Ht and we denote it simply: u < on dO. More generally, if M is a random 
field defined on [0, T] x O, we note u < M on dO if u — M < on dO. 

From now on, we can follow step by step the proof of the maximum principle for OSPDE 
in the homogeneous case in |llj : the first step consists in establishing an estimate for the 
positive part of the solution with null Dirichlet condition. To get this estimate, we can 
adapt to our case the arguments of proof of Proposition 5.2 in then Ito formula for 
the difference of 2 elements in TZi oc (^, /, g,h, S) (Proposition 5.3 in This yields the 

comparison theorem (see Theorem 5.4 in 

Theorem 6. Assume that dO is Lipschitz. Let (£*, /*, g, h, S l ), i = 1,2, satisfy assumptions 
(H), (HIL), (OL) and (HOL). Consider («*,!/*) G n Xoc (f , /*, g, h, S 1 ) , i = 1,2 and 
suppose that the process (n 1 — u 2 ) + belongs to TIt an d that one has 

E^\\f^.,., U 2 ,Vu 2 )-f^.,., U 2 ,V U 2 )\\* # . t y<oo, for all t€[0,T\. 

If t 1 < S, 2 a.s., f 1 (t,uj,u 2 ,Vu 2 ) < f 2 (t,uj,u 2 ,Vu 2 ), dt <g> dx <g> dP-a.e. and S 1 < S 2 , 
dt ® dx <S> dP-a.s., then one has v}{t, x) < u 2 (t, x), dt <S> dx ® dP-a.e. 



By adapting the proof of Theorem 5.5 in we get first the maximum principle in the 
case u < on dO: 

Theorem 7. Assume that dO is Lipschitz and suppose that Assumptions (H), (OL), 
(HOL),(HI2p), (HOoop) and (HD6lp) hold for some 9 G [0, 1[, p > 2 and that the 
constants of the Lipschitz conditions satisfy 

B 2 

a + Y + 72 P 2 < A - 
Let (u, v) G TZioc (^, /, g, h, S) be such that u + G 7i. Then one has 



£ KL.oo.-t * k ^ E ( - S °\L + (\\f°'% t r + (HlsTQ 5 + (P°l a U 5 

+ ii(^)iiL+(||/'i;/+(iib'i 2 ii;^+(iii^i 2 ii;^) 

where k (t) is constant that depends on the structure constants and t G [0,T]. 

As in the homogeneous case (see Theorem 5.6 in we can generalize the previous result 
by considering a real Ito process of the form 

oo ft 



ft ft 

M t = m+ b s ds + y2 <Tj,sdBi 
Jo =1 Jo 



r- 

where m is a random variable and b = (bt)t>o, o = (oi,t, cr nit , ...)t>o are adapted pro- 
cesses. 
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Theorem 8. Assume that dO is Lipschitz and suppose that Assumptions (H), (OL), 
(HOL),(HI2p), (HOoop) and (HD(9p) hold for some 9 G [0,1[, p > 2 and that the 
constants of the Lipschitz conditions satisfy 

B 2 

a + y + 723 2 < A. 

Assume also that m and the processes b and a satisfy the following integrability conditions 



E \m\ p < oo, E 



\b s ds 



p(l-0) 



p(i-fl) 



< oo, E [ / \a s \ i-s ds 



< oo, 



for each t G [0, T\. Let (u,v) € 7Zi oc (£, f, g,h, S) be such that (u — M) + belongs to TLt- 
Then one has 



< c(p)k(t)E[ - m )+ - (S' - m)||^ + 

+ \\(S'o -mm 



+ 



£ 

* \ 2 



=0|2 



+ (\\(f'-b) + 



+ 



e 

* \ 2 



:o 1 2 



Hiu^i 2 y 2 + ii^--i 2 fc v 



(16) 



where k (t) is the constant from the preceding corollary. The right hand side of this es- 
timate is dominated by the following quantity which is expressed directly in terms of the 
characteristics of the process M, 



c(p)k(t)E[\\(C - m)+ - (S' - m)||^ + {\\f°'% t ) + 



=0|2 



+ 



h 



0|2 



E 

* \ 2 



+ ll^o -m) + \\l + 



+ 



/' 

Jo 



\bJ i-e ds 



p(i-0) 



+ 



+ 



|cr.,|i-9 ds 



p 

* \ 2 



+ 



p(l-B) 



\h'\ 



p 

* \ 2 
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